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We propose a spin manipulation technique based entirely on electric fields applied to acceptor states 
in p-type semiconductors with spin-orbit coupling. While interesting in its own right, the technique 
can also be used to implement fault-resilient holonomic quantum computing. We explicitly compute 
adiabatic transformation matrix (holonomy) of the degenerate states and comment on the feasibility 
of the scheme as an experimental technique. 
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The physical realization of quantum computing rests 
on the ability to reversibly manipulate two level systems 
called qubits. While the promise of high computational 
power is certainly a tantalizing one, the intrinsic chal- 
lenges associated with decoherence, adiabatic evolution, 
control and noise errors in quantum gate operations are 
still to be mastered. 

One ingenious way to overcome quantum noise errors is 
the use of Non-abelian Holonomic (Geometric) Quantum 
Computation schemes Q. In these procedures, through 
the slow tuning of some external parameters such as ap- 
plied magnetic or electric fields, the qubit evolves adi- 
abatically (with constant energy) around a path that 
changes its eigenstate from an initial to a final state. 
Generically, this quantum evolution is free of dynami- 
cal factors and is geometric in nature, depending only 
on the path in parameter space. Geometric holonomy 
could constitutes a fault-tolerant way to perform quan- 
tum computation Q, Although some experimental 
systems that would exhibit such behavior have been pro- 
posed 0, , holonomic quantum computing overall still 
lacks the variety of concrete application proposals that 
conventional quantum computation enjoys. 

In this paper we propose using electric fields to manip- 
ulate the spin of acceptor states in semiconductors with 
spin-orbit coupling such as Ge, GaAs and Si. The ac- 
ceptor impurity ion will bind a p-hole from the spin 3/2 
valence band of the semiconductor [|| and the full hamil- 
tonian of the impurity system in an electric field is given 
by the linear or quadratic Stark effect Q. Spin-orbit 
coupling is essential to the existence of the Stark effect 
in these semiconductors. There are two doubly degener- 
ate Kramer states for any value of the electric field and 
slowly rotating the electric field induces SU (2) rotations 
in the degenerate eigenstates of each energy level. When 
the electric field swaps over a cycle and returns to its ini- 
tial orientation the holonomy matrix is dependent on the 
geometry of the swap only. Consistent with prior the- 
oretical analysis, we take these holonomies to represent 
quantum gate transformations 

We begin by a short introduction of the main idea 
of holonomic quantum computing. We then introduce 
the hamiltonian of the acceptor states in an electric field 
and show how any SU (2) holonomy can be obtained by 



changing the field's orientation adiabatically, hence pro- 
viding the basis for a set of gate operations. We shortly 
discuss the viability of such a scheme and the experimen- 
tal challenges involved. We close by proposing an alter- 
native scheme, using external uniaxial strain, which can 
be used to achieve spin manipulation in the absence of 
any external field. The holonomic qubits discussed here 
are based on the same principle as the recently discov- 
ered dissipationless spin current in hole doped semicon- 
ductors |2(- The individual qubits can therefore be cou- 
pled to each other by a quantum bus architecture based 
on the dissipationless spin current, offering exciting new 
possibilities towards the realization of an all solid-state 
holonomic quantum computer. 

In the thesis of holonomic quantum computing, quan- 
tum information is encoded in an n-fold degenerate 
Hilbert space of a hamiltonian H\ dependent on some 
external 'control' parameters (fields) A |||. Upon a 
cyclical change of these parameters around a loop C 
during time T such that A, n = \ ou t, the system will 
evolve between the initial state (n vector) \ip)i n into 
\i>)out = e le ° T r(C)|'0)i„, where eo is the initial eigenvalue 
H\ in \ip)i n = eo\tp)i n . The first factor is just the dynami- 
cal phase, and will be omitted, while the second factor is 
the non-abelian Wilczeck-Zee jjjj curvature connection 
(matrix) : 

r(C)=Pe^o^ 5 A £ =(^(A)|^-|^'(A)>, (1) 

where i, j — 1, ...,n and where P represents the path or- 
dering due to the fact that the gauge connection A 1 - 1 is 
now a matrix acting on the degenerate space of Hamilto- 
nian eigenstates (p, denotes the different control param- 
eters). The degenerate Hilbert space of the Hamiltonian 
encodes the quantum information where the eigenstates 
are the codewords while the non-trivial holonomies asso- 
ciated with it represent the unitary transformations or 
'computations' over the code. Zanardi and Rasetti 
showed that this prescription is sufficient to implement 
quantum computation on single qubit holonomic gates. 
In subsequent papers [1 01 [HI U3 several schemes 
for realizing holonomic computation have been proposed. 
The schemes involve geometric manipulation of trapped 
ions |l2j . charge pumping within Josephson junction net- 
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works 0, and Josephson charge qubits 0. Controlled 
manipulation of U(l) holonomies (Berry phases) using 
nuclear magnetic resonance on a system of weakly cou- 
pled H 1 and C 13 nuclei has been experimentally achieved 
with great accuracy by Jones et al [l3|. 

It would be of great advantage to have a conventional 
solid-state system where holonomic computation can be 
implemented by using only electric fields. In this pa- 
per we are concerned only with single qubit holonomy, 
leaving multi-qubit ones for a later publication 0] ■ We 
look at p-type cubic symmetry semiconductors such as 
Ge, GaAs and Si. The strong spin-orbit coupling present 
in these systems breaks up the valence bands into two 
doubly degenerate bands of spin 3/2 with helicity ±1/2 
and ±3/2. The double degeneracy is nothing else than 
Kramers degeneracy and is guaranteed by T-invariance. 
T-invariance is maintained even when acceptor impuri- 
ties (B, Al, Ga, In) are introduced in the semiconductors. 
The holes that bind to these impurities will maintain a 
certain symmetry subgroup of the original cubic symme- 
try of the valence bands they came from. Let us now 
consider the effect of an applied external electric field 
E on the acceptor-bound hole state. For large electric 
fields, the field distortion near the impurity ion can be 
safely neglected and the acceptor-hole state has the cubic 
symmetry of the crystal Td x /, giving rise to a quadratic 
Stark effect 0: 
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where £j is the ionization energy, p — er the dipole mo- 
ment (r being the mean radius of the ground state), and 
S are the spin- 3/2 matrices, describing the valence band 
states, which essentially have P 3 / 2 character . We have 
also defined {A, B} = (AB + BA)/2. Readers familiar 
with semiconductor theory will recognize in the form of 
H E 2 the Luttinger Hamiltonian structure, with the sub- 
stitution k — > E. This is no coincidence since the sym- 
metry group of both Hamiltonians is the same. For small 
applied electric field, we must take into consideration the 
local field of the ions, thereby reducing the symmetry 
from Td x / to Td and giving rise to a linear Stark effect 
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-lj={E X {Sy, S Z }- 



■E y {S z ,S x } + E z {S x ,S y }), (3) 



where p = ea E with a E the Bohr radius. The constants 
a,/3, 5, x,r are given in Table 1, although the estimates 
for x m the literature vary considerably (x — 0.26 ac- 
cording to Kopf and Lassman so the value in Table 
1 should be taken as a lower limit). We want to mention 



TABLE I: Values of the coefficients in the linear and quadratic 
Stark Hamiltonians 





a 


P 


5 


X 


r(A) 


Ge 


1 


-0.3 


-0.36 


0.7 x 10 -3 


91 


Si 


1 


-0.2 


-0.42 


1 x!0~ 2 


34.4 



that the donor and acceptor Hamiltonians and physics 
are essentially different, with the donors undergoing only 
a quadratic Stark shift as opposed to the acceptor com- 
bination of the above linear and quadratic shifts. 
Although for some field E the acceptor Hamiltonian will 
be a weighted sum of linear and quadratic Stark effects, 
we prefer, without any loss of generality, to work in cither 
of the two regimes and not in the intermediate one. Each 
of the Hamiltonians above has two doubly degenerate 
Hilbert spaces, roughly corresponding to values of the z- 
component of the spin S z being either ±1/2 or ±3/2 (this 
would be exactly true if the Hamiltonians were isotropic) . 
The 'control' parameters are the components of the elec- 
tric field E. We must now show we can achieve 'quan- 
tum computations over the code'. These are represented 
by SU (2) holonomies over each degenerate Hilbert space 
(equivalently, we must show that we can move within an 
energy subspace by adiabatically changing E). 

We now prove that such holonomies do indeed exist 
in our system and give an explicit generic procedure to 
calculate them. While we could just brute-force diag- 
onalize the Hamiltonians above and treat each of them 
separately, we prefer to use a more elegant approach that 
reveals more of the Hilbert space structure. This was de- 
veloped by Demler and Zhang 0] in the context of the 
SO (5) theory of high Tc superconductivity, and extended 
by Murakami, Nagaosa and Zhang to the case of hole 
band in semi-conductors. Readers not interested in the 
derivation can jump to the next page where we give the 
expression for the SU(2) holonomies. Out of the spin-3/2 
J x , J y , J x we can define the new 4x4 matrices: 



n/3 



{S y ,s z }, r 



{s z ,s x }, r 3 



n/3 



{Sy, S x } 



r 4 = -^-s*), r 5 = sf-^/ 4x4 , (4) 

which satisfy the SO(5) Clifford algebra r<T b + r fc r a = 
2<5 a ;,i4x4. In explicit form, these matrices are: 




where a 1 , i — 1,2,3 are the usual Pauli matrices and / 
is the identity matrix ( 2 x 2 in this case). We observe 
that the Hamiltonians He and H E 2 can now be cast into 
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a new clean form 



H E = d E I + d a E T a - H E 2 = d% 2 I + d a E2 T a ; a = 1, 5 

(6) 

where / is the 4x4 identity matrix and where the fol- 
lowing identities hold: 

d E = d E = d% = 0, d% = pxE x , d% = p X E y , d% = PxE z 



2 2 2 

d° E 2 = aE 2 , d E 2 = 5E z E y , d E 2 = SE Z E X 

Sj, Ej_ Ej_ 



d E 2 — ——SE x E y , d E 2 



d%,=- P fb(2El-El-El) 
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Since the two Hamiltonians for small and large field have 
been brought to the same symbolic form, we can manip- 
ulate them together and only substitute for the values of 
d a at the end of the calculation. The deep physical reason 
as to why the two apparently different Hamiltonians are 
actually very similar is the unbroken T-invariance of the 
system that leads to Kramers degeneracy. The eigenval- 
ues are e± = d°±d where d — \J d a d a and they depend on 
the electric field. For the linear Stark effect explicit sub- 
stitution shows us that the split is independent on the 
direction of the electric field E while for the quadratic 
case, it strongly depends on its orientation. The Clifford 
matrices T a have two eigenvalues, each two-fold degener- 
ate (this is obvious from the form of T 5 ). The gauge con- 
nection then represents SU (2) adiabatic changes on the 
two-fold degenerate sub-bands ± hence the total gauge 
group is SU(2)^ x SU(2) + = SO(4). Identical to the 
work of Murakami Nagaosa and Zhang ^jj) we can define 
the projection operators into the two energy subspaces 
H = e+P+ + e-P~: 

p+ = i ( i + ^n, p- = ld- d in- (8) 

Adiabatic rotation of the field E implies moving within 
one of the subspaces of energy e. We can define a covari- 
ant gauge field strength: 

r)P+ r)P~ 1 

where T ab = ^[r a ,r b ] are the generators of the 50(5) 
algebra. The field A a lives in the space of the d a, s but 
our control parameters are the electric field components 
Ei. We hence have to do a 'coordinate transformation' 
and obtain: 



which gives us a holonomy computation when the electric 
field is varied between Einmai = Efi„ a i over an arbitrary 
(closed) curve C 

r(C) = P exp( j> A l dE t ) = P exp(- I ^J^^^) 

\ (11) 

Let us, without any loss of generality momentarily fo- 
cus on the e+ subspace. By choosing specific rotations 
(specific contours C) of the field Ei we can change an 
initial state \tp)i n — (1,0,0,0) into the degenerate state 
within the same energy level, i.e. \ip}out — r(C)|V')m — 
(0, 1,0,0). In fact, in the general case, starting from an 
arbitrary \tp)i n we can reach, through carefully choosing 
the contour C, any other eigenstate within the degenerate 
subspace by electric field manipulation. In a physical in- 
tuitive picture, the spin within the e + subspace will follow 
the electric field as it tries to stay within the energy sub- 
space. We have hence achieved spin manipulation with 
electric fields and showed that holonomic computation is 
possible in semiconductors with spin-orbit coupling. 

In general, due to the non-abelian nature of Ai, the 
path ordered integral has to be done numerically, over in- 
finitesimal segments in parameter space and taking into 
account that different components of A do not commute 
with each other. While this is more of a nuisance than an 
intellectual challenge, it is comforting to know that for 
certain curves the expression can be simplified and path 
ordering can be easily implemented while still maintain- 
ing the full capability to transform the eigenstates into 
one another. We give such examples for both the linear 
and the quadratic Stark effect below. 

For the linear Stark effect, again working in the e + 
energy subspace the expression for the holonomy Te(C) 
becomes particularly simple: 



T E {C) =Pexp(- 
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2d 2 dEi 
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where the are the 3 pauli matrices. In polar coor- 
dinates E — (E sin8 coscj), -E sin (9 sin Ecos9) for con- 
tours C which keep constant the absolute value of the 
electric field, we find that spherical triangles between the 
points A (0 = 0,0 = 0i), B ( 6 = tt/2,0 = fa), and C 
(6 = 7r/2,0 = fa) are particularly easy to path order. 
Since we are changing only one angle at a time achieving 
this technologically should be easier than trying to imple- 
ment variations in both angles (although, as Zee points 
out, there is a bit of confusion on how to go 'around the 
corners' |18j 1 

For the case of the quadratic Stark effect things are 
more complicated. While finding a nice form for the 
holonomy factors in the general case is almost impos- 
sible due to the anisotropy in the Hamiltonian He^, we 
can look at the idealized spherical symmetric situation 
for which [3 = 8/y3. This does not introduce large er- 
rors, as the anisotropy in these materials, although sig- 
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nificant, is still small enough so that the spherical ap- 
proximation works well. In this case we find, in units of 
-f|: d a E2 T a = [3{E ■ S) 2 - 5/4pE 2 h xi . The holonomy 
structure resides exclusively in the first term. In fact, 
with the electric field replaced by a magnetic field, this is 
exactly the Hamiltonian studied by Ze e 1181 in explain- 
ing a pioneering experiment by Tycko [2l|. The gauge 
field in polar coordinates is = cos 0cr 3 /2 — sin do\ and 
Ae = (7 2 . For spherical triangles starting at 9 = go- 
ing to some value 9 on an arc of fixed <f>\ (which we can 
choose to be zero for convenience) then going at fixed 
9 on an arc to some non-zero <f> and then back to the 
north pole along constant longitude, the holonomy reads 
I>(C) = W^VW [Hi where: 

W{~ 1 = exp(—i9(cos4>a2 — sin^cri)), W = exp(z<72#) 



V — exp(— i— cr 3 ) exp(i — (cos0CT3 — 2sin6>cri)) (13) 

We now turn to the problem of the feasibility of the 
scheme proposed for spin-manipulation by the Stark ef- 
fect. We need that the coherence time of spins of bound 
holes be larger than the time in which we can adiabat- 
ically rotate the electric field. New experiments showed 
that the coherence time is larger than 1ms 0, [2oj |. 
justifying the use of acceptor-bound-hole wavefunctions 
as qubits. It is indeed difficult to perform experiments 
which probe non-abelian phase factors. The original work 
of Tycko |2l| and subsequently the more complete exper- 
iment by Zwanziger, Koenig and Pines on nuclear 
magnetic quadrupole resonance proved the existence of 
the Wilczeck and Zee non-abelian transport of degener- 
ate states. Instead of rotating the applied fields and keep- 
ing the sample fixed, these experiments kept the applied 
field fixed and rotated the sample, which is an equivalent 
procedure. The rotation frequencies were of the order 
of a few kHz (2020 Hz in Zwanziger et al). Imagining 
an electric-field version of this experiment, the rotation 
period of the field is already less than the bound hole 
spin coherence time, but further improvement may be 
necessary for a realistic measurement. We also need to 
guarantee, during the field rotation, that the adiabatic 
approximation is accurately maintained and that the ac- 
ceptor impurity is not ionized. The ionization energies 
for acceptor states are of the order 10 — 60 meV (see 
Table 2). The splitting between the two levels e+ and 
e_ can be computed from our expressions for their en- 
ergies (using the constants in Table 1) and are of the 
order 10 meV for Ge in a field of 10 6 V/m. Hence the 
ionization and splitting energies are roughly the same 
size and much larger than the applied electric field frc- 
quecy of rotation. The frequencies required for varying 
the electric field are hence low enough as to cause nei- 
ther ionization of the impurity-hole system nor a break- 
down of adiabaticity. The variation of dynamical phases 



TABLE II: Ionization energies in meV for different acceptor 
impurities (B, Al, Ga) in Si and Ge 





B 


Al 


Ga 


Ge 


10.4 


10.2 


10.8 


Si 


45.0 


57.0 


65.0 



over the sample volume, which usually leads to extensive 
dephasing can be overcome by an electric field variant 
of the double-sweep spin-echo techniques which refocus 
inhomogeneities in the dynamical phase but double the 
effect of the geometric phase . 

In conclusion, this paper presents a novel way to ma- 
nipulate the spins of acceptor impurity-bound hole states 
in p-type semiconductors with strong spin-orbit coupling 
using electric fields. Depending on its magnitude, the 
electric field couples both linearly and quadratically to 
the spin of the acceptor state through the Stark effect, 
but although apparently opposite, the two effects still 
maintain the T-invariance of the underlying semiconduc- 
tor. The spin manipulation is completely geometric and 
realizes, in a practical solid-state system, the theoretical 
proposal for holonomic quantum computing. We have 
obtained an explicit and general form for the holonomy 
matrix which transforms adiabatically transports degen- 
erate eigenstates. While our analysis is specific to spin 
3/2 it is trivially generalized for any spin, provided the 
Stark effect is present. We have also briefly analyzed 
the experimental feasibility of the scheme. In a future 
work[l^, we explore the idea of different qubits commu- 
nicating via a bus architecture based on the dissipation- 
less spin current 

In closing, we want to stress that the wide variety of 
options available for spin manipulation in semiconductors 
with spin-orbit coupling makes them particularly suitable 
for the realization of quantum computers. Aside from 
electric field manipulation, the spin of the acceptor states 
(as well as that of the free holes) is heavily influenced by 
interior and exterior applied strain on the material, so 
that the same analysis presented above can be done for 
strain-induced holonomies and spin-manipulation. 

In the closing stages of this work, we noticed the 
independent, recent work by Yuri Serebrennikov which 
presents similar ideas to the ones exposed here [23J. 
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